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Abstract 

A review of interesting current topics in turbulence 
research is decorated with examples of popular fallacies 
about the behaviour of turbulence. Topics include the 
status of the Law of the Wall, especially in compressible 
flow; analogies between the effects of Reynolds number, 
pressure gradient, unsteadiness and roughness change; the 
status of Kolmogorov’s universal equilibrium theory and 
local isotropy of the small eddies; turbulence modelling, 
with reference to universality, pressure-strain modelling 
and the dissipation equation; and chaos. Fallacies include 
the mixing-length concept; the effect of pressure gradient 
on Reynolds shear stress; the separability of time and space 
derivatives; models of the dissipation equation; and chaos. 

1. Introduction 

I first met Keith Stewartson in the early 1960’s, when 
he was a young member of the British Aeronautical Re- 
search Council’s Fluid Motion Subcommittee and I was its 
(very) young secretary. Even in those days, I dimly sensed 
that Keith was not particularly fond of turbulence. It is, 
therefore, a matter of double regret that I should be giv- 
ing, so soon, a lecture in his memory, and should be forced 
to choose the subject of turbulence as being my only area 
of aerodynamic competence. 

Those who knew Keith will recall that his strongest 
term of scientific condemnation was “unrigorous” . Pm sure 
he regarded the whole phenomenon of turbulence as being 
unrigorous and probably invented by the Devil on the sev- 
enth day of Creation (when the Good Lord wasn’t look- 
ing); I am inclined to agree. Keith would certainly have 
approved of the rigour of Horace Lamb’s “Hydrodynam- 
ics” (Cambridge University Press) - what the reviewer of 
a later book once called his “awful correctness”. Lamb, 
after discussing all the branches of hydrodynamics known 
to him, finally had to deal with turbulence and remarked, 
in Article 365, p. 651 of the 1916 edition, “It remains 
to call attention to the chief outstanding difficulty of our 
subject.” Seventy-odd years have come and gone; difficul- 
ties in hydrodynamics have come and gone; but turbulence 
still remains as the “chief outstanding difficulty of our sub- 
ject”. Another dead friend, Jack Nielsen, Chief Scientist of 
NASA Ames, said a few years ago that turbulence model- 
ing was the “pacing item” in the use of the NAS computer 
complex, and I think his comment, like Lamb’s, is still 
true. 

In the last ten years or so we have become able to solve 
the complete time-dependent Navier-Stokes equations for 
turbulent flow. However, the Reynolds numbers at which 


usually only about three or four times the lowest at which 
turbulence can exist, and are considerably lower than the 
Reynolds numbers obtainable in laboratory experiments, 
let alone those found in real life. Therefore, although 
turbulence is starting to become accessible to computers, 
there is no immediate prospect of the subject going the 
same way as stress analysis and succumbing almost en- 
tirely to computation: unlike elasticity, turbulence is a 
non-linear (strictly, quasi-linear) phenomenon and, at least 
at high Reynolds numbers, is at present accessible only to 
experiment. Thus, experimental fluid dynamics will last 
for many years (hopefully, for my working lifetime). 

Of course, turbulence would merely be a laboratory 
curiosity or a computational playground if it were not for 
its extreme importance in real life and in all the scien- 
tific and engineering disciplines represented here today 
in meteorology, aeronautical aerodynamics, shipbuilding, 
oceanography, in all forms of pipeline design and manufac- 
ture, in combustion, in any form of mixing of contaminant, 
whether of heat or concentration or pollutant - in other 
words in almost all forms of “interesting” fluid motion ex- 
cept those on an extremely small scale. The cream poured 
into a cup of coffee goes turbulent, and the flow patterns 
look very cloud-like, (The poem on the letter “H” will be 
quoted in the oral lecture.) 

I propose to use this Memorial Lecture to try to inject 
a certain amount of rigour into the study of turbulence, 
specifically by using the occasion to review some popular 
fallacies about turbulence and the way in which turbu- 
lent flows behave. Some of these fallacies or illogicalities 
are propagated by popular but outdated textbooks, but 
some are at a deeper level of incomprehension, including 
the preconceptions of workers in statistical mechanics who 
think that turbulence must be easy. Naturally, parts of 
the material that I will produce are controversial, in the 
sense that some of my professional colleagues may disagree 
with me. However, I hope that even the controversial sec- 
tions of the paper will be of interest and may stimulate 
clarifying discussion, either at this meeting or after it. It 
is of course difficult to group illogicalities into any logical 
order, so I have imbedded them into a study of the more 
popular topics of turbulence “theory”. I hope the result is 
neither a rag-bag nor a grab-bag. The oral lecture will be 
less specialised than this written version. 

My favorite definition of turbulence is that it is the 
general solution of the Navier-Stokes equations. This is 
the perfect answer by a government servant to an inquiry 
by a Congressman or Member of Parliament: it is brief, it 
is entirely true, and it adds nothing to what was known 


already. Nearly everybody believes, of course, that the 
Navier-Stokes equations are an adequately exact descrip- 
tion of turbulence, or indeed of any other nonrelativis- 
tic motion of a Newtonian fluid. Even the smallest ed- 
dies in turbulence in ordinary liquids and gases at earth- 
bound temperatures and pressures are large compared to 
the mean free path between molecular collisions, so the 
constitutive equation of the fluid is not in doubt. How- 
ever, Sec. 4 of the present paper deals with the influence 
of fluctuating dilatation div u on turbulence in compress- 
ible gas flow, and in this case the uncertain value of the 
bulk viscosity 0 (Goldstein 1 ) may matter. 

Fortunately for professional educators, it is generally 
accepted that the basic phenomena of turbulence are the 
same at any Mach number - except for some special effects 
to be discussed in Sec. 4 - so unless stated otherwise I will 
assume the density to be constant. 

2. The Law of the Wall 


One of the main building blocks, or even foundation 
stones, of the engineering study of turbulence is the “Law 
of the Wall”. It derives from the hypothesis / assumption 
that, sufficiently close to a solid wail (meaning, for exam- 
ple, a distance from the wall an order of magnitude less 
than the diameter of a pipe or the thickness of a boundary 
layer) the flow depends only on the distance from the wall, 
on the shear stress at the wall t w , and on fluid properties. 
The characteristics of the outer part of the flow do not 
matter except that they determine r w . (In the discussion 
below, the term “shear stress” will sometimes be used to 
mean “shear stress / density”, for short.) 

Let us consider a boundary layer for simplicity. The 
characteristics of the outer part of the flow to be consid- 
ered include the free-stream velocity U e and the boundary 
layer thickness 6. The irrelevance of U ei as such, is a conse- 
quence of Galilean (translational) invariance and does not 
need much discussion, The irrelevance of 6 is more crucial, 
as it depends on the assumption that the flow close to the 
surface consists of eddies whose length scales (in all direc- 
tions) are proportional to y, with negligible contributions 
from eddies whose length scales (in any direction) depend 
on 6: if this is so, the boundary layer thickness should not 
appear in any scaling of the inner-layer eddies. We shall 
see in Sec. 5 that this hopeful view is not quite correct, 
but it is certainly acceptable to first order. 

The consequence of these arguments is, of course, that 
the mean velocity and turbulence near the surface should 
scale on the “friction velocity” u r = (r w /p) 1 / 2 , on the 
distance from the surface, y, and on the kinematic viscosity 
v. One of the several dimensionally-correct ways of writing 
this relationship is 

U/u T = fi{u T y/u). (1) 

Another \b obtained by differentiating Eq. (l) and hiding 
a factor of u r yjv inside the function / 2 , as 


dU/dy = (u r /y)/ 2 (u r y/t/). (2) 

Here u r yju is an eddy Reynolds number based on the eddy 
velocity and length scales, i.e. the friction velocity and the 
distance from the surface. At large values of this Reynolds 
number we expect the effects of viscosity on the turbulence 
to be negligible and therefore Eq. (2) reduces to 

dU/dy = u T /(*y) (3) 

where k « 0.41 is a constant - Von Kar man’s constant, of 
course. The integral of this relationship is the logarithmic 
law, the additive constant C « 5 being a constant of inte- 
gration depending on the velocity difference between the 
wall and the point at which Eq. (3) becomes valid. 

The advantages of the above analysis over the tradi- 
tional “overlap” demonstration are (i) that the only as- 
sumption made about the outer layer is that it doesn’t 
matter, and (ii) that a simple physical argument can be 
used to simplify Eq. (2) to the so-called mixing-length 
formula, Eq. (3). 

The constant of integration C is equal to 5 only on 
smooth walls: on rough walls, it becomes a function of 
the roughness Reynolds number xi r k/ v and of the rough- 
ness geometry; the uncertainty of the effective origin of 
y on rough walls is a further complication. The constant 
k, on the other hand, is supposedly universal: it is the 
same in flows of water and of air on all geometries involv- 
ing smooth surfaces, and indeed on all geometries involv- 
ing only small roughness; it is the same in the atmospheric 
boundary layer, in the depths of the ocean and on the sands 
of Mars. Alas k and C are not constant within the tur- 
bulence modelling community - a remarkably wide range 
of values is in use. Those quoted we from the painstaking 
data analysis of Coles 2 . 

Now there are still textbooks - and even living people 
- that regard the log law as a deduction from the mixing- 
length formula, Eq. (3), (which it is) and also regard 
the mixing-length formula for the inner layer as correct 
(which it is) and also regard Prandtl’s original derivation 
of the mixing-length formula by analogy with molecular 
motion as correct (which it certainly is not) . As the Roman 
Catholic Church quite properly pointed out to Galileo, the 
success of deductions from a hypothesis does not prove 
its truth. Philosophers call this the fallacy post hoc , ergo 
propter hoc (“after that, therefore because of that”) and 
it is the basis of witch-doctoring (last time we slaughtered 
a white cow, it rained; there is a drought; therefore. . .). 

Quite apart from philosophical questions of falsifiabil- 
ity, it is clear that if a result can be derived by dimen- 
sional analysis alone, like Eq. (3), then it can be derived 
by almost any theory, right or wrong, which is dimen- 
sionally correct and uses the right variables. There is a 
strong suspicion that Prandtf got the idea of the lumps of 
fluid (“Flussigkeitsballen”) of mixing-length theory from 
visual studies of turbulent open-channel flows with parti- 
cles sprinkled on the surface to show up the motion. Unfor- 
tunately the boundary condition at a free surface permits 


only motion tangential to the surface and not normal to it, 
so the surface becomes a plane of symmetry with the vor- 
ticity vector everywhere normal to it. The only motions 
that can remain are what sailors, but not landlubberly tur- 
bulence researchers, call “eddies”. Try it, and you will see 
what Prandtl saw. 

3. Extensions to the law of the wall 

The law of the wall derived in Sec. 2 is valid, or is 
supposed to be valid, for a shear stress equal to the wall 
shear stress and a density equal to the wall density. There 
is some support for an extended version of Eq. (3), still for 
u r ylu > 30 approx., in conditions where either the shear 
stress r = —puv or the density p varies with distance from 
the surface. If u r is replaced by (r/p) 1 / 2 , we get 

dUfdy = (r//») 1/a /(«v) ( 4 ) 

The hand-waving argument for Eq. (4) is that, in the orig- 
inal analysis leading to Eq. (3), u r is really being used as 
the scale at height y, and not as a true surface parameter: 
if r varies with y then the local value, rather than the wall 
value, is the correct one to use in formulating an eddy ve- 
locity scale. This would be a rigorous argument only if the 
typical eddy size were small compared with y, so that the 
local shear stress would be closely equal to the right basis 
for a velocity scale, namely some kind of weighted-average 
shear stress over a y distance equal to a typical eddy size. 
Unfortunately, of course, the eddy size is of the same order 
as y. 

All we can claim is that local shear stress gives the 
best easily-available velocity scale. Therefore, the exten- 
sion of Eq. (3) to Eq. (4) requires an extension of faith in 
the inner-layer hypothesis which by no means all research 
workers possess. Nevertheless the application of Eq. (4) to 
flows with suction or injection, where the shear stress varies 
with distance from the surface according to r = t w + pUV w 
is quite well supported by experiment. An operational dif- 
ficulty is that in typical flows with suction or injection the 
surface is porous, on a length scale h, say, which is usually 
not small compared with the viscous scale i//u T , so that 
the “roughness” or “porosity” Reynolds number u T hju is 
important, implying that the additive constant in any in- 
tegral of Eq. (4) will depend on the surface conditions as 
well as on the transpiration parameter V w /u r . 

4. Compressible flow 

In the inner layer of a boundary layer in compressible 

flow, the shear stress is approximately equal to the surface 

value, but the density varies quite rapidly with distance 

from the surface (increasing as the temperature decreases 

with distance from the hot wall). The “Van Driest trans- 

formation” transforms inner-layer velocity profiles to fit 
the incompressible log. law. The transformation is, in 
effect, an integral of Eq. (4) with p as a function of y. 
Here T and hence p come from the assumption of a con- 
stant turbulent Prandtl number: details will not be given 
here, but can be found in Ref. 3 and elsewhere. The Van 
Driest skin-friction formula is derived from the Van Driest 
transformation. Predictions of skin friction in compressible 


boundary layers (on flat plates in zero pressure gradient, 
say) are currently a subject of controversy, but there are 
certainly no experimental data that reliably invalidate the 
Van Driest skin-friction formula or the Van Driest trans- 
formation. This is probably the best justification for the 
extension of the law-of- the- wall analysis discussed in Sec. 
3, but doubtless does little for the confidence of the deter- 
minedly subsonic. 

In low-speed flow, the mean (streamwise) pressure 
gradient, as such, has almost no effect on turbulence (see 
Sec. 6). In compressible flow, streamwise pressure gradi- 
ents change the density of fluid elements and can produce 
large changes in turbulence quantities, especially, of course, 
in flows through shock waves (e.g. Selig et al. 4 * * ). Moreover, 
even pressure fluctuations which are not small compared 
with the meaji pressure can affect turbulence. Specifically, 
if the Mach number based upon a typical fluctuating ve- 
locity and the local speed of sound is no longer small com- 
pared to unity, there may be significant dissipation of tur- 
bulent energy via dilatation fluctuations dttm, and signif- 
icant correlations between fluctuations of pressure and of 
dilatation 5,6 . Measurements correlated by Birch k Eggers 7 
show that the rate of spread of a turbulent mixing layer 
(in zero mean pressure gradient) starts to depend signifi- 
cantly on Mach number at Mach numbers close to unity. 
The more more recent data of Papamoschou k Roshko 8 * 
show even larger Mach-number dependence. This appar- 
ently contradicts the well-known finding that the behav- 
ior of compressible boundary layers can be quite well pre- 
dicted by turbulence models that ignore compressibility 
effects (except of course that the right mean density must 
be used), at least for Mach numbers up to about 5. How- 
ever, the typical turbulence intensity of a mixing layer is 
about five times that in a boundary layer, which implies 
that a mixing layer at M=l, where M is based on the 
mean velocity difference across the layer, has the same ra- 
tio of velocity fluctuation to speed of sound (a.k.a. fluctu- 
ating Mach number) as a boundary layer at roughly M=5. 
There is great current interest, stemming from the NASP 
and SCRAMJET projects, in prediction of mixing layers 
as the only shock-free turbulent flow for which the data 
show obvious effects of compressibility. 

5. “Inactive” motion 

The log-law analysis relies on the first-order hypothe- 
sis that u r , y and v are the only relevant variables, which 
cannot be exactly and perfectly true. If the arguments 
that lead to Eq. (3) are applied to the turbulent motion 
they lead to results for the log-law region like t ?/u 2 = 
constant, whereas any boundary-layer experiment shows a 
decrease with increasing y, starting as close to the wall as 
u r y/i/ = 17 at typical small laboratory Reynolds numbers. 
This has led some people to regard the whole law-of-the- 
wall concept of local scaling as fallacious and its apparent 
success for the mean motion as fortuitous. Fortunately, 
this apparent discrepancy in the log-law analysis can be 
used to rescue the basic assumptions, by taking note of 
the so-called “inactive” motion ®’ 10 . The concept is sim- 
ple: the motion near the surface, even though it results 
mainly from eddies actually generated near the surface, 


is necessarily affected by eddies in the outer part of the 
flow (i.e. those whose length scale is of the order of £), 
Because the pressure fluctuation at a given point in a tur- 
bulent flow is derived from an integral of the governing 
Poisson equation over the whole of the flow, it follows that 
the eddies in the outer part of the boundary layer or pipe 
flow can produce pressure fluctuations which extend to- 
wards the surface and cause nominally-irrotational motion 
in the surface layer. An equivalent, alternative, explana- 
tion is the “splat” mechanism (the origin of the term will 
be explained in the oral lecture) in which the large eddies 
in the outer flow are supposed to move towards the surface, 
to be reduced to rest by the normal-component “imperme- 
ability” condition at the wall, and to release their normal- 
component energy into the two tangential components u 
and w. 

The “splat effect” motions, and the pressure fluctu- 
ations generated in the outer layer, have very long wave- 
lengths in the x and z directions compared to the motions 
generated close to the surface. It follows from the con- 
tinuity equation that the v-component velocity produced 
near the surface by outer-layer pressure fluctuations or 
large-eddy intrusions is of the order of y/X times the u- or 
u;-component velocity, where A is the x- or z-component 
wavelength. Therefore the contribution of the “inactive” 
motion to the shear stress —put; is small, of the order of 
y/X - hence the name “inactive”. Note that the “inac- 
tive” fluctuations are not entirely irrotational: the bound- 
ary condition u = 0, v — 0 at the surface results in the 
generation of a Stokes layer (see Sec. 6 on “slip veloc- 
ity”). Even though “inactive” u-component fluctuations 
contribute significantly to u 2 , producing the anomalous in- 
dependence of u 2 mentioned above, the effect on the mean 
law of the wall is very small. (A logarithm is a slowly 
changing function, so that fluctuations in u r have very lit- 
tle effect on the term ln(u T y/i/) in the log. law, and, there- 
fore, the time-average velocity closely follows the log. law 
written with time-average u T .) The same arguments can 
be used to support the use of the log. law in unsteady-flow 
calculations at not-too-high amplitudes. The unsteady log. 
law must also be limited to not-too-high frequencies of un- 
steadiness: one would expect it to break down, at given y, 
at a frequency which was not small compared to the typical 
turbulence frequency u T /y. Very few unsteady-flow exper- 
iments reach frequencies high enough to disturb the log. 
law - which is a criticism of unsteady-flow experiments in 
general. 

The contribution of the “inactive” fluctuations to the 
power spectra of u and w at low wave numbers (low fre- 
quencies: wave number = 2tt/ [wavelength]) is consider- 
able, resulting in very large differences between the mea- 
sured spectra in typical turbulent flows and those predicted 
by inner-layer analysis. The latter predicts that the wave- 
number spectral density should scale on u r and y, and that 
the wave number k should appear as Icy (since we have ne- 
glected i/, this applies only for u r y/v > 30 and at wave 
numbers small compared with the viscous limit, but nei- 
ther restriction concerns us here). In practice, there is an 
apparent Reynolds-numb er effect at given u T y/i/: strictly 
it is a y/6 effect, but y/6 = [u r y j u) j {u r S j u , 


In the atmospheric boundary layer, which is of the or- 
der of 1 km thick, the inactive-motion effects on spectra 
measured at the standard height of 10 m are very large, and 
in particular the u-component spectrum follows a —5/3 
power law down to very low wave numbers. This phe- 
nomenon, which is present, but less spectacular, in lab- 
oratory boundary layers, has been the cause of a large 
amount of confusion, controversy and difficulty, because 
the classical Kolmogorov scaling indicates that the spec- 
trum should vary as only for wave numbers large 

compared to those of the energy-containing eddies. In the 
context of the atmospheric boundary layer at a height of 10 
m this means wavelengths much smaJler than 10 m. The 
fact that the experimentally-observed spectrum follows the 
—5/3 law down to wave numbers far lower than could be 
expected from the arguments of inner-layer scaling and the 
Kolmogorov universal-equilibrium hypothesis is one of the 
most difficult “fallacies” in turbulent flow: it is of course 
a case of post hoc ergo propter hoc . 

In s umm a r y, the qualitative idea of “inactive motion” 
explains both the apparent failure of inner-layer scaling 
and the unexpected success of the -5/3 law. 

6. “Slip velocity” 

Several difficulties or misconceptions about turbulent 
flows over walls can be cleared up if we recall that the 
very thin viscous wall region u r y/y < 30 really produces 
what might be called a “slip velocity” between the fully- 
turbulent flow and the surface. As well as the obvious 
example of Reynolds-number (and Pec let-number) effects, 
they include the effects of pressure gradient, unsteadiness 
and change of surface roughness. 

6.1 Effects of Reynolds number and Peclet number 
(viscosity and conductivity) 

If the Reynolds number of a turbulent flow - based 
on total thickness and, say, the square root of the maxi- 
mum shear stress or turbulent energy - is large, classical 
(e.g. Kolmogorov) theory suggests that the details of the 
turbulent motion should be independent of Reynolds num- 
ber, except for the very smallest eddies which are responsi- 
ble for viscous dissipation of turbulent kinetic energy into 
thermal internal energy. In this respect at least, classi- 
cal theory seems to be correct, and there is no significant 
evidence to refute it. If the Reynolds number of a given 
turbulent eddy, made with its typical velocity fluctuation 
and its typical length scale, is large, there is no reason why 
viscous effects on the eddy should be significant. (This 
statement should strictly be phrased in statistical terms!) 
In a pipe flow, half the mean-square u-component inten- 
sity near the centre-line comes from wavelengths larger 
than the pipe diameter, so the “eddy Reynolds number” 
of the main energy-containing eddies is of the same or- 
der as the mean Reynolds number defined at the start of 
the paragraph, and we can use the former for simplicity. 
The “energy cascade” process of Kolmogorov theory, at- 
tributable to random vortex stretching, implies that tur- 
bulent energy is transferred from energetic eddies of low 
wave number (i.e. large Reynolds number) to weak ed- 
dies of high wave number (small Reynolds number), and 



although back-scatter transfer from small eddies to large 
can occur intermittently, the time-average transfer of en- 
ergy is from the large eddies to the small and there seems 
to be no significant “back scatter” of viscous effects. 

Near a solid surface (y + > 30) the largest eddies, 
whose wavelength is roughly equal to y, are no longer very 
large compared to the smallest eddies (the smallest-eddy 
scale, Kolmogorov’s rj or /*, is about 0.06y at y + = 30), so 
the energy-containing eddies - which also carry the shear 
stress - start to depend on viscosity. (Also, and slightly 
differently, the mean velocity gradient becomes so large 
that viscous shear stress is a significant fraction of the to- 
tal shear stress.) Therefore, viewed from the outer part of 
the flow, there is a viscosity-dependent region near the wall 
and so the velocity difference between the surface and, say, 
y/S = 0.1, depends on Reynolds number. Viewed from 
the outer part of the flow, there is a Reynolds-number- 
dependent “slip velocity” at (strictly near) the surface. 

In a free shear layer (wake, jet, mixing layer. . .) there 
is no true viscous effect unless the Reynolds number is 
so low that turbulence can only just exist. However, free 
shear layers can be quite strongly dependent on the ini- 
tial conditions, for long distances downstream, and since 
the initial conditions frequently do depend on Reynolds 
number there is a “pseudo-viscous” effect. 

A corollary of the negligibility of viscosity as part of 
turbulent transport of momentum is the negligibility of 
conductivity in the transport of heat or mass by turbu- 
lence. Briefly (again) the “turbulent Prandtl number” is 
independent of the molecular Prandtl number unless the 
Reynolds number based on eddy velocity scale and eddy 
length scale, i.e. u r y/i/ is small. 

6.2 Effect of pressure gradient 

Another of the standard incomprehensions about tur- 
bulent flow is the effect of (streamwise) mean pressure gra- 
dient on the turbulence as such: recall that we are con- 
sidering only incompressible flow. It arises partly because 
experimenters tend to normalize their turbulence measure- 
ments by the local mean velocity. In adverse pressure gra- 
dient, say, the mean velocity decreases with increasing x 
so the normalized turbulence intensities, shear stress etc. 
increase, However it can easily be shown that absolute tur- 
bulence properties on a given streamline are only slightly 
affected by pressure gradient. 

The Reynolds-s tress transport equations do not con- 
tain the mean pressure (they contain correlations between 
the pressure fluctuation and instantaneous rate of strain, 
but pressure fluctuations have no connection whatsoever 
with the mean pressure). Also, the x-component mean 
vorticity dVjdx - dU/dy is unaffected by pressure gradi- 
ent, and if we assume that the boundary layer approxi- 
mation is valid this means that dU/dy is unaffected, even 
though the pressure change leads to a change in veloc- 
ity all through the shear layer and thus may change S 
significantly. Alternatively, recall that a static- pressure 
gradient does not affect the total pressure P (on a given 
streamline) directly. In the simple case of two-dimensional 


flow, therefore, dP/dx/j, where 0 is the stream function, 
is unaffected, and if dp/dy is negligible, as required by 
the boundary-layer approximation, a little algebra shows 
that dU/dy is unaffected. At the surface, where the total 
pressure is equal to the static pressure, there is a change 
in P and dU/dy , produced of course by viscous stresses. 
The “internal layer” , in which the total pressure and mean 
vorticity rise to their unaffected profiles, gradually spreads 
out from the surface, but outside this the static-pressure 
gradient has no effect except to reduce the mean velocity 
and thus thicken the boundary layer. This result applies to 
laminar or turbulent boundary layers (or other wall flows 
such as those in tapered ducts). In summary, the initial 
effect of pressure gradient is confined to the “slip velocity” 
at the wall. 

Mean pressure gradients do have some effects on the 
turbulent motion. Adverse pressure gradient stretches ed- 
dies in the y direction, because the shear layer thickens: 
however, the area, in side view, of a given eddy or fluid 
element is unaltered, and so if we suppose that the length 
scale of an eddy is just the square root of its area in side 
view, or the cube root of its volume, the length scale is 
unaltered. (This is admittedly a crude argument.) Of the 
terms in the Reynolds-stress transport equations, the only 
ones directly affected are the y-component diffusion terms, 
which are the derivatives of various triple prodects, etc., 
with respect to y. If the triple product on a given stream- 
line is unaffected but the streamlines diverge in the x - y 
plane because of the adverse pressure gradient, the y-wise 
derivative is reduced. 

6.3 Unsteadiness 

The effect of unsteadiness can be understood in the 
same way as that of pressure gradient - of course, unsteadi- 
ness is usually forced by a streamwise pressure gradient. 
In the case of unsteady laminar flow the internal layer is 
called a Stokes layer. There are close correspondences in 
lammar flow between an infinite oscillating plate in still 
air and flow over an infinite stationary surface driven by 
an oscillating pressure gradient, and the qualitative corre- 
spondence carries over to turbulent flow. If the pressure 
gradient is strong enough to cause separation (however de- 
fined), the internal layer is carried into the outer part of 
the flow and the “slip velocity” concept breaks down, as it 
would in steady separation. 

6.4 Change of roughness 

Another occasion where a change of boundary con- 
dition affects the flow only in an “internal layer” is the 
flow downstream of a change in surface roughness. This 
is comparatively rare in aerodynamics but an important 
case in meteorology where, for example, air can flow from 
the “smooth” ocean to the land and undergo a change 
of apparent surface roughness. Indeed, the internal- layer 
concept was first proposed to describe this case. As the 
surface boundary condition changes, the additive constant 
C in the logarithmic law for a smooth surface is replaced 
by the appropriate value for a rough surface. The effect of 
this change in surface boundary condition spreads outward 
from the surface at an angle of the order of r ms u/27, i.e. 


of the order of 3 %, so that the rate of contamination of 
outer-layer turbulence by inner-layer changes is no greater 
than about 1 or 2 degrees. Since the pressure gradient is 
nominally zero there is no streamline divergence above the 
internal layer, although the change in velocity in the inter- 
nal layer produces a vertical displacement of the outer flow 
(upwards, in the case of a smooth-to-rough change where 
the flow In the internal layer is retarded). 

7, Spectra and convection velocity 

Classical turbulence theory aims to predict all the sta- 
tistical properties, not simply the Reynolds stresses. In 
particular it deals with the statistical distribution of eddy 
sizes. It is usually formulated in terms of wave-number 
spectra, wave number being a vector with the direction of 
wavelength and the magnitude of 2tt / wavelength. (The 
alternative is two-point spatial correlations, which are less 
convenient mathematically.) Wave-number spectra are the 
Fourier transforms of the two-point correlations, but a full 
description requires correlations for all magnitudes and di- 
rections of the distance between the two points, or spectra 
for all magnitudes and directions of the wave number). In 
most experiments only frequency spectra, and a few cor- 
relations along the coordinate axes, are measured. 

This is the best place to comment on the definition 
of “frequency” in turbulence. The frequency seen by an 
observer moving with the mean flow is (velocity scale of 
turbulence) / fJength scale of turbulence ) - for example 
u T /y in the inner layer - but the frequency seen by a fixed 
observer is approximately (MEAN velocity) / (length scale 
of turbulence) and is usually much larger. The recipro- 
cal of the moving-observer frequency is sometimes called 
the “eddy turnover time”: this is of course an order-of- 
magnitude concept. A related difficulty is the status of 
time derivatives: all transport equations in fluid flow, in- 
cluding the Navier-Stokes equations, have the operator 


on the left-hand side. It is called the substantial deriva- 
tive, or the transport operator, and it is the rate of change 
with time seen by a fluid element. The relative size of the 
temporal and spatial derivatives depends on the velocity 
of the observer but the sum of the derivatives does not. 

The fixed-observer frequency is used to deduce x-com- 
ponent wave-number spectra from frequency spectra, using 
Taylor’s hypothesis that the speed at which the turbulence 
pattern moves downstream (its “convection velocity”) is 
closely equal to the mean velocity. It is qualitatively ob- 
vious that this will only work well if the mean velocity is 
large compared to the velocity scale of turbulence, so that 
an eddy is carried past the measurement point in a time 
very much less than its turnover time. A more precise 
analysis is possible. 

There are various definitions of the actual “convec- 
tion velocity” of turbulence: most are in effect phase ve- 
locities and therefore not ideal for considering convection 
of turbulent kinetic energy or Reynolds stress. A plausi- 


ble definition of a group (energy-transport) velocity comes 
from considering the streamwise (say, x-wise) “diffusion” 
of turbulent energy (transport of the turbulent energy by 
the turbulence): the energy flux rate, whose x derivative 
app ea rs in the turbulent energy equation, is p*u/p + (u 3 + 
uv 2 + uw 2 )/ 2. Rates of energy flux due to pressure fluc- 
tuations seem to be small - except perhaps near the free- 
stream edge of a turbulent flow where pressure fluctuations 
drive an “irrotational” motion which extends outside the 
vortical region - and are certainly not measurable, which 
is some justification for neglecting them. Doing this, and 
writing q 2 for u 2 + v 3 + u/ 2 , (so that the turbulent kinetic 
ene rgy is q 2 / 2) , the above energy flux rate can be written 
as (u 3 -j- uv 2 + uw 2 )/q 2 . We can define the transport ve- 
locity of turbulent energy as this flux rate divided by the 
turbulent energy. The largest contribution to the numera- 
tor is « 3 /2 - though the others are not negligible - so the 
transport velocity is of order \/(q 2 ) x S u , where S u is the 
skewness of u. Now S u lies in the range ±1 approx, over 
most of a boundary layer, so we can finally say that the 
x-component transport velocity of turbulent energy is not 
more than a few times >/(q 2 )' Since this is the difference 
between the group velocity of the turbulence and the mean 
velocity, we see that the difference is a small percentage of 
the mean velocity in flows with low turbulence intensity, 
such as boundary layers. This quantitatively justifies the 
use of Taylor’s hypothesis in such flows and of course allows 
an estimate of its inaccuracy in highly- turbulent flows. 

Differences between convection velocity and mean ve- 
locity are large near the free-stream edges of mixing layers 
and jets. In these regions the irrotational motion, induced 
by pressure fluctuations generated in the high-intensity re- 
gion of the flow near the inflexion point (s) in the velocity 
profile, is strong compared to the true (vorticity-carrying) 
turbulence, and its convection velocity is necessarily close 
to the mean velocity in the high-intensity region. The ro- 
tational motion (vorticity pattern) seems to travel at a 
speed close to the local mean velocity, as predicted by the 
above analysis (intensities near the outer edge of a jet are 
not large). In terms of the above analysis, the streamwise 
transport velocity of the vorticity pattern is still domi- 
nated by the triple-product terms, while ~p r u/p determines 
the transport velocity of irrotational motion. 

The de Havilland Comet I jet airliner had four engines, 
buried in the wing roots. The designers carefully arranged 
that the jets themselves would clear the fuselage, but forgot 
the “near field” pressure fluctuations - far more intense 
than the jet noise - that drive the irrotational motion. 
The pressure patterns, travelling at the above-mentioned 
convection velocity, produced fluctuating stresses at the 
fixed-observer frequency in the aircraft skin, which led to 
fatigue of the aluminium. 

Later marks of Comet had the engines toed out. 

Misconceptions about turbulence can be expensive! 

8. The microscale and the Kolmogorov theory 

Frequently, the Taylor “microscale” is used as a length 
scale in discussions of wave-number (or frequency) spectra. 


The microscale A is a hybrid scale of turbulence. It is 
usually defined by 

A 2 = u 2 /(du/dx ) 2 (6) 

(other definitions with different choices of velocity compo- 
nent or gradient direction occasionally appear). This is 
an equation whose numerator is a property of the energy- 
containing turbulence, but whose denominator is a prop- 
erty of the dissipating eddies (if the dissipating eddies ar e 
statistically isotropic the dissipation rate is lSu(du/dx) 2 ). 
For this reason it is a misconception to regard the mi- 
croscale as the length scale of any particular group of ed- 
dies: it actually lies closer to the length scale of the dis- 
sipating eddies than that of the energy-containing eddies. 
The Reynolds number based on the microscale and the 
root-mean-square turbulence intensity,)^) 1 / 2 /!/, however, 
has a more understandable meaning. If the Reynolds num- 
ber is high enough for the dissipation to be equated to 
the isotropic formula, the microscale Reynolds number is 
proportional to the square root of an “eddy* Reynolds 
number for the energy-containing motion, based on the 
rms turbulence ..intensity and the dissipation length scale 
t = ( u 2 ) 3 / 2 /£. Of course, this does not give the microscale 
the status of Eddy Length Scale post hoc . 

It is important to notice that the “dissipation" in the 
definition of L is in fact the rate of transfer of turbulent 
kinetic energy from the large eddies to the smallest eddies 
which is, by all prevailing turbulence theories, supposed to 
be a property of the large eddies rather than the smallest 
eddies. The smallest eddies simply rearrange themselves 
to dissipate the energy handed down to them. If the turbu- 
lence is changing slowly with time (or streamwise distance) 
then, of course, the rate of transfer from the large eddies 
to the smallest eddies is equal to the rate at which en- 
ergy is being dissipated by the smallest eddies, but this 
is not formally an equality because the “cascade” process 
is not instantaneous. In rapidly-changing turbulent flows 
the “equilibrium” arguments fail, and the rate of transfer 
from the energy-containing eddies to the dissipating eddies 
is not equal to the rate at which energy is being transferred 
from the dissipating eddies to heat. 

This restriction on Kolmogorov’s “universal equilib- 
rium” theory, which we used in Sec. 6.1, is too often for- 
gotten. 

Another restriction of the Kolmogorov theory is that, 
of course, energy which is transported in the y direction 
by turbulent “diffusion” will be generated at small y, but 
dissipated at large y where the statistical properties are 
different. In particular, in flows with a free-stream bound- 
ary, energy is generated in regions of large mean shear and 
then transported in the positive y direction to regions of 
zero or negligible mean shear before being dissipated. The 
energy transfer through the inertial subrange at the second 
location is likely to be intermediate between the dissipation 
rates at the two locations. 

Nevertheless results from a large number of experi- 
ments on turbulent shear layers have recently been anal- 
ysed 11 to show that Kolmogorov scaling works remark- 


ably well when adjusted for the intermittency factor 7 (the 
fraction of time for which the flow at a given location is 
turbulent). In an intermittent region, the average of any 
turbulence quantity within the turbulent part of the flow 
is 1/7 times the conventional average over all time. For 
example the conventional-average spectral density and the 
dissipation c must both be multiplied by 1/7. However 
the Kolmogorov “-5/3” law for the spectral density in 
the so-called inertial subrange contains e 2 / 3 so that, for- 
mally, there is a spare factor of 7 1 / 3 and we certainly do 
not expect the Kolmogorov law to hold if written with 
conventional-average quantities. The data analysis of Ref. 
11 shows that the Kolmogorov formula still works for a 
wide range of intermittent flows when written for the tur- 
bulent part of the flow, i.e. taking account of the “spare 
factor”, and using the dissipation rate at the local value 
of y. Since the formula strictly applies only to nearly- 
homogeneous turbulence, and an intermittent region, al- 
most by definition, contains only one large eddy at a time, 
this result is a surprising testimonial to the robustness of 
the Kolmogorov theory. Needless to say, the usual cautions 
about post hoc apply. 

9. Turbulence modelling 

9.1 Normal pressure gradients 

An incomprehension entirely unrelated to turbulence, 
which nevertheless causes confusion in tests of turbulence 
models, is the effect of normal pressure gradient on bound- 
ary layers and other shear layers. If the shear layer obeys 
the boundary layer approximation then, by definition, the 
pressure gradient in the y direction is negligibly small. 
However, if in a real flow the normal pressure gradient 
is not negligible, there will be a velocity gradient dU/dy 
even in the external stream (where the total pressure is 
constant) and this velocity gradient will, in principle, lead 
to extra production of turbulence via the product of mean 
velocity gradient and turbulent shear stress. Of course, the 
same effects would be found within the shear layer y < 6 , 
but would be less easily identified. Therefore, even if a 
turbulence model produces exactly correct predictions of 
the shear stress - given the mean velocity profile as input 
- it will not give acceptable results in the case where nor- 
mal pressure gradients affect the mean velocity gradient. 
(Recall that the boundary- layer momentum equation can 
be written as dP/dx = dr/dy.) This is probably a much 
more important reason for inaccuracy of predictions based 
on the boundary layer approximation in rapidly-growing 
flows near separation than the often-quoted presence of 
significant normal-stress gradients. 

9.2 Universality 

Perhaps the biggest fallacy about turbulence is that 
it can be reliably described (statistically) by a system of 
equations which is far easier to solve than the full time- 
dependent three-dimensional Navier-Stokes equations. Of 
course the question is what is meant by “reliably”, and 
even if one makes generous estimates of required engineer- 
ing accuracy and requires predictions only of the Reynolds 
stresses, the likelihood is that a simplified model of tur- 


bulence will be significantly less accurate, or significantly 
less widely applicable, than the Navier-Stokes equations 
themselves - i.e. it will not be “universal”. 

Irrespective of the use to which a model will be put, 
lack of universality may interfere with the calibration of 
a model. For example, it is customary to fix one of the 
coefficients in the model dissipation-transport equation so 
that the model reproduces the decay of grid turbulence 
accurately. This involves the assumption that the model is 
valid in grid turbulence as well as in the flows for which it 
is intended - presumably shear layers, which have a very 
different structure from grid turbulence. 

It is becoming more and more probable that really 
reliable turbulence models are likely to be so long in devel- 
opment that large-eddy simulations (from which, of course, 
all required statistics can be derived) will arrive at their 
maturity first. (The late Stan Corrsin once described the 
process of turbulence modelling as a “trek to determi- 
nacy”.) Certainly, over the last twenty years the rate of 
progress in turbulence modelling has been pretty small 
compared to the rate of progress in development of dig- 
ital computers, and the consequent increase in Reynolds- 
number range and geometrical complexity attainable by 
simulations. Until recently, most work has concentrated 
on “complete” simulations, covering the whole range of 
eddy sizes, while large-eddy simulations, which alone of- 
fer the prospect of predictions at high Reynolds numbers, 
have been somewhat neglected. 

9.3 Eddy viscosity and gradient transport 

Turbulence models which invoke an eddy viscosity (of 
whatever type) necessarily produce pseudo-laminar solu- 
tions with the stresses closely linked to the mean-flow gra- 
dients: they may be well-behaved but they are not usually 
very accurate away from the flows for which they have 
been calibrated. Turbulence models based on term-by- 
term modelling of the Reynolds-stress transport equations 
produce solutions which may be accurate in some cases, 
but are liable to fail rather badly in other cases: that is, 
they are “ill-behaved” in a way that eddy-viscosity meth- 
ods are not. 

It may be this “reliable inaccuracy”, rather than the 
larger computer resources needed for Reynolds-stress trans- 
port models, which has led to two-equation (e.g. jfc,e) 
or even one-equation methods being the industry stan- 
dard. With all goodwill to my friends Barrett Baldwin 
and Harv. Lomax, the one-equation Baldwin-Lomax tur- 
bulence model has been extended - by others - far beyond 
its intended domain, simply because it has the virtue of 
almost never breaking down computationally! 

It has, of course, often been said that it is just as un- 
reliable and unrealistic to define an eddy viscosity entirely 
in terms of turbulence properties (as in the k, e method) as 
to define it entirely in terms of mean-flow properties as in 
the Baldwin-Lomax method. Eddy viscosity is the ratio of 
a turbulence quantity (i.e. a Reynolds-stress) to a mean- 
flow quantity (i.e. a rate of strain or velocity gradient), so, 
like the microscale, it is a hybrid quantity. 


Minor fallacies in turbulence modelling abound, but 
misuse of gradient-transport hypotheses is probably re- 
sponsible for more than its fair share. One of the most 
spectacular was the use many years ago, by authors I will 
not identify, of the gradient-transport approximation for 
diffusion of turbulent energy by pressure fluctuations. In 
terms of classical physics, anything less likely than pres- 
sure diffusion to obey a gradient-transport approximation 
could scarcely be imagined. A fallacy which has, in charity, 
to be regarded as a deliberate approximation, is the use - 
even in Reynolds-stress transport models - of the eddy- 
diffusivity (gradient-transport) approximation for the tur- 
bulent transport terms. It appears that most of the tur- 
bulent transport of Reynolds stress is provided by triple 
products of velocity fluctuations, rather than by the pres- 
sure diffusion just mentioned, and therefore a gradient- 
transport approximation is not so obviously unphysical. 

9.4 The dissipation-transport equation 

Most turbulence models, whether relying on an eddy 
viscosity or on the Reynolds-stress transport equations, 
use the dissipation-transport equation to provide a length 
scale or time scale of the turbulent flow. Strictly, the 
length scale or time scale required is that of the energy- 
containing Reynolds-stress-bearing eddies, not that associ- 
ated with the dissipating eddies as such, and so two ques- 
tions arise. One is whether the rate of dissipation is ade- 
quately equal to the rate of energy transfer from the large 
eddies (which clearly, is the quantity that we really want 
to model); the other is whether, if we really pretend to 
be using the dissipation transport equation - all of whose 
terms depend on the statistics of the smaiiest eddies - 
we can logically model those terms by using the scales 
of the larger, energy-containing eddies. I think it is in- 
escapable that current models of the so-called dissipation 
transport equation, which certainly do parameterize the 
terms as functions of the large-eddy scales, start out with 
the dissipation-transport equation as such and end up with 
a totally-empirical transport equation for the energy trans- 
fer rate. In other words, the relation between the “dissipa- 
tion” transport models and the exact transport equation 
for turbulent energy dissipation is so tenuous as not to 
need consideration. Unfortunately, even Reynolds-stress 
transport models usually employ this suspect dissipation- 
transport equation to provide a length scale, and this is 
undoubtedly one of the reasons why Reynolds-stress trans- 
port models have not outstripped two-equation models. A 
less-used alternative to the £ equation is the w equation 
(admitted to be totally empirical), u is nominally pro- 
portional to t/k where k is the turbulent kinetic energy, 
but conversion from one to the other (in either direction) 
produces the interesting result that the turbulent trans- 
port terms in the transport equation for the first quantity 
(the integral of transport terms over the flow volume be- 
ing by definition zero) convert to a transport term plus 
a “source” term in the equation for the second quantity. 
There is increasing evidence that using w to provide a 
length scale gives better results than using t: if there is 
a reason other than more judicious choices of empirical co- 
efficients, it must lie in the above-mentioned source term. 



9.5 Invariance 


One of the customary requirements of a turbulence 
model is that it should be “invariant” (with respect to 
translation or rotation of axes). The boundary layer (thin- 
shear layer) equations are not invariant: it is therefore 
quite unrealistic to expect a shear-layer model to be totally 
invariant, and it is perfectly realistic to suppose that the 
direction normal to the shear- layer (y) is a special direc- 
tion. There seems to be no reason why a turbulence model 
should not, given an identifiable “special direction” in a 
shear- layer use that special direction for orientation of its 
empirical constants and functions. Even though equations 
(such as the Navier-Stokes equations or the time-average 
Reynolds equations) may be invariant, the boundary con- 
ditions for which they are to be satisfied certainly are not 
invariant (almost by definition). Therefore, the solutions 
of the exact, or approximate, equations of motion of turbu- 
lent flow cannot be expected to be invariant with respect to 
translation or rotation. From this it is a rather small step 
to argue that the empirical constants or functions in these 
model equations should, again, be released from invariance 
requirements. 


9.6 Local modelling of pressure-fluctuation terms 

The mean products of fluctuating pressure and fluc- 
tuating rates of strain that act as redistribution terms 
in the Reynolds-stress transport equations represent, very 
crudely speaking, the effect of eddy collisions in making 
the principal Reynolds stresses more nearly equal - that 
is, making the turbulence more nearly isotropic (statisti- 
cally). The shear stress in isotropic turbulence is zero, so 
the effect of the pressure-strain terms on the shear stress, 
and their modelling, is of great interest. 


Pressure fluctuations within a turbulent flow are one 
of the Great Unmeasurables: they are of the order of 
pu 2 and so, unfortunately, are the pressure fluctuations 
induced on a static-pressure probe by the velocity field. 
That is, the signal-to-noise ratio is of the order of one. To 
say that signals cannot be educed even with S/N = 0(1) 
is itself a fallacy, but in this case the attempts made to do 
so have not met with general acceptance. Pressure fluc- 
tuations can be extracted from simulations, but these are 
confined to low Reynolds number. 


An equation for the pressure (mean and fluctuating) 
can be obtained by taking the divergence of the Navier- 
Stokes equations. It is a Poisson equation, and it is nec- 
essary in turbulence modelling to consider the different 
terms on the right-hand side separately, by writing a Pois- 
son equation for each and adding the solutions to get the 
pressure. One such is the equation for the “rapid” pres- 
sure, which for a two-dimensional boundary- layer flow is 


VV 2 dUdv 

p dy dx’ 


( 7 ) 


The “rapid” pressure is so called because it responds im- 
mediately to a change in the mean flow, as represented by 
dU/dy. To regard this apparently-surprising fact as phys- 
ically meaningful is a misconception: it is just a result of 


the way we take averages, and, obviously, the turbulence 
at a given instant does not know what the mean flow is. 
A highly symbolic solution of the equation is 


£=-2V- 2 ^ 
p dy dx 


( 8 ) 


where V 2 is a weighted integral over the whole flow vol- 
ume. In other words, the “rapid” pressure at a given point, 
and its contribution to the pressure-strain terms at that 
point, depend on conditions for a distance of several typical 
eddy length scales around that point - i.e. they are “non- 
local”. The same non-locality accounts for the presence 
of irrotational velocity fluctuations outside the turbulent 
motion. 


Almost all current stress-transport turbulence models, 
with the exception of that of Durbin 12 , model the pressure- 
strain terms and other pressure- velocity correlations en- 
tirely as functions of local quantities. (All the other terms 
in the Reynolds-stress transport equations are genuinely 
local quantities.) This is equivalent to replacing Eq. (8) 
by 


V — = 

P dy dx 


( 9 ) 


that is, evaluating dU/dy at the position where p 7 is 
required and volume-integrating only dv/dx. 


In Ref. 13, the behavior of existing models for the 
pressure-strain terms was analyzed, using simulation data 
in a duct flow to evaluate the terms directly. The results, 
surprisingly, suggest that the difference between the ex- 
act pressure-strain terms, using p' from Eq. (8), and the 
approximate results, using p ,B from Eq. (9), is negligibly 
small (or, at least, small enough to be hidden in the empir- 
ical coefficient in the pressure-strain model) except in the 
viscous wall region. Within the viscous wall region, the 
difference between the true pressure fluctuation p 7 and the 
approximate pressure fluctuation p ,B is not only very large 
but eccentrically behaved. It is not suggested that viscous 
effects, arising from the v = 0 boundary condition at the 
surface, are directly to blame: it is much more likely that 
the effects of the v = 0 boundary condition are mainly re- 
sponsible, but it is surprising that these effects should be 
small outside the viscous wall region. A final possibility 
is that the changes in turbulence structure with t i T yjv in 
the viscous wail region are so large as to invalidate local 
models. 


This suggests not only that standard pressure-strain 
models are grossly inaccurate in the viscous wall region, 
but also that any extension of a standard turbulence model 
into the viscous wall region will be similarly inaccurate. 
This inaccuracy can be camouflaged by the insertion of 
“low-Reynolds-number" functions, nominally functions of 
the wall distance u T yj i/. Obviously, if the real flow scales 
with u r y/i/, this simple procedure suffices, but if the flow 
approaches, or goes beyond, separation then inner-layer 
scaling - and presumably “low-Reynolds-number” models 
- break down, even if u T y/u is replaced by the guaranteed- 
real quantity k l / 2 y/v. 


10. Chaos 

“What kept you?” you may ask. Chaos has been 
one of the buzzwords in applied mathematics in recent 
years, and turbulence is often cited as the supreme ex- 
ample. The complication of turbulent motion, with its 
broad spectrum of wavelengths, is far greater than that of 
the “chaotic” solutions of some low-order systems of cou- 
pled ordinary differential equations. Analysis of simulation 
data 14 suggests that the dimension of the turbulence at- 
tractor (roughly, the number of modes or “degrees of free- 
dom” needed to represent the turbulent motion) is several 
hundreds at least, even at the lowest Reynolds number at 
which turbulence can exist. The upper bound on the di- 
mension is, roughly, the number of totally-arbitrary modes 
(say, Fourier modes or finite-difference formulae) needed to 
represent the motion. Now since direct-simulation calcu- 
lations need, typically, 128 3 » 2 x 10 6 Fourier or finite- 
difference points for flows at a very modest laboratory- 
scale Reynolds number, we can take the upper bound of 
the attractor dimension as being of this order: for the 
barely- turbulent flow of Ref. 13, 32 s « 30000 might do. 
Large-eddy simulations need fewer points: 128 3 might do 
for any Reynolds number, at least if the viscous wail region 
did not have to be resolved. These are all impracticably 
large estimates of the attractor dimension. 

However, several authors have based their work on the 
classically incorrect syllogism “Solutions of some equations 
with few degrees of freedom yield complicated behavior: 
turbulence has complicated behavior: therefore turbulence 
may be represented by the solution of equations with few 
degrees of freedom” . The last hypothesis of course stood 
by itself for many years B.C. (before chaos), and a great 
deal of brain power has been applied to prove it - i.e. to 
produce a usably small set of modes to describe turbulence 
- but without great success: the most ambitious efforts 
require an amount of computing time which is not much 
less than that of a large-eddy simulation. 

The concepts of chaos theory may of course be quali- 
tatively useful in turbulence studies. One is the concept of 
predictability. Qualitative arguments about the non-linear 
Navier-Stokes equations suggest that if two almost identi- 
cal turbulence fields with the same boundary conditions 
are set up at time t — 0, then the two instantaneous veloc- 
ity and pressure fields will become more and more different 
at time goes on, even though the statistical properties of 
the two fields will still be (nearly) equal. To a worker 
in turbulence, particularly an experimenter, this does not 
seem odd - but the issue of instantaneous versus statisti- 
cal predictability has attracted a lot of attention in chaos 
studies, and perhaps our intuition about the Navier-Stokes 
equations may be put on a firmer footing. Deissler 15 re- 
views applications of chaos studies in fluid dynamics; for 
a popular introduction to chaos studies in general, see the 
book by Gleick 16 ; and see also, of course, the new inter- 
disciplinary journal “Chaos’. 

11. Conclusions 

In this paper, we have gone all the way from very basic 


questions of turbulence theory to the important practical 
question of the reliability of turbulence models, and then 
ended in chaos. The fallacies that we have discussed do not 
necessarily form a coherent story, but I think it can be said 
that most of them fall into the general category of wishful 
thinking - the hope of finding simple solutions to a difficult 
problem. I will end with one of my favorite quotations, 
from H. L. Mencken, “to every difficult question there is a 
simple answer - which is wrong” . 
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It is to be noticed that the ratio of (8) to (12) U of the order y/(ac/t>) t numerical factors 
being omitted. In all cases to which our approximations apply this ratio is large, so that 
the radial vibrations are much more slowly ertinguiahed, so far as viscosity alone is con- 
cerned, than those which correspond to values of n greater than 0. This is readily accounted 
for. In the latter modes the condition that there is to be no slipping of the fluid in contact 
with the vessel implies a relatively greater amount of distortion of the fluid elements, and 
consequent dissipation of energy, in the superficial layers of the gas. 

The method of the dissipation function, which was applied in Art. 348 to the case of 
water waves, might be used to obtain the result (12) for the radial vibrations, but would 
lead to an erroneous result for « >0, since the underlying assumption that the motion is 
only slightly modified by the friction is violated at the boundary. 

In the gravest radial vibration we have ka ■ 4-403, whence 

t = *0743-. 

V 

In the case of air at 0° C. this makes r * -fiOa 1 *. 


Turbulent Motion. 

365. It remains to call attention to the chief outstanding difficulty of 
our subject. 

It has already been pointed out that the neglect of the terms of the 
second order (u5u/3z, & c.) seriously limits the application of many of the 
preceding results to fluids possessed of ordinary degrees of mobility. Unless 
the velocities, or the linear dimensions involved, be very small the actual 
motion in such cases, so far as it admits of being observed, is found to be 
very different from that represented by our formulae. For example, when 
a solid of 'easy’ shape moves through a liquid, an irregular eddying motion 
is produced in & layer of the fluid next to the solid, and a trail of eddies is 
left behind, whilst the motion at a distance laterally is comparatively smooth 
and uniform. 

The mathematical disability above pointed out does not apply to cases 
of rectilinear flow, such as have been discussed in Arts. 330, 331;" but even 
here observation shews that the types of motion investigated, though always 
theoretically possible, become under certain conditions practicallv unstable. 

The case of flow through a pipe of circular section was made the subject 
of a careful experimental study by Reynolds f, by means of filaments of 
coloured fluid introduced into the stream So long as the mean velocity 
( w o) over the croBssection falls below a certain limit depending on the radius 
of the pipe and the nature of the fluid, the flow is smooth and in accordance 

• This Art. is derived with alight alteration from a paper cited on p. 636. 

t “An Experimental Investigation of the Circumstance* which determine whether the 
Motion of Water shall be Direct or Sinuous, and of the Law of Resistance in Parallel Channels.” 
Phil. Trans, t. cIxxjt. p. 935 (1883) {Papers, t. il p. 51]. For a historical account of the researches 
and partial anticipations of other writers, see Knibbs. Proc. Roy. Sot. R’.S.W. t. rxxi. p, 314 
(1897). Reference is made in particular to Hagen, BerL Abh. 1854, p. 17. 



